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ABSTRACT 

The radiative decays D° — > K°P 0/ -f with nonresonant K°P° (P° = tt°, rj, rj') are con- 
sidered within the framework which combines heavy quark effective theory and the chiral 
Lagrangian. Due to neutral mesons the amplitudes do not have bremsstrahlung contri- 
butions. We assume factorization for the weak matrix elements. Light (virtual) vector 
mesons are found to give the main contribution to the decay amplitudes. The decay 
D° — > K°TT 0/ -f is predicted to have a branching ratio of 3 x 10~ 4 , with comparable contri- 
butions from parity-conserving and parity-violating parts of the amplitude. The decays 
with 77(77') in the final state are expected with branching ratios of 1.1 x 10~ 5 and 0.4 x 1CT 7 
respectively and are mainly parity conserving. 



Radiative weak decays of mesons have proved to be a fertile ground for the investigation 
of models of the strong interactions involved, and of the nonleptonic weak Hamiltonian. 
This is true for long-distance dominated K decays |]J as well as for B radiative decays, 
driven by short-distance dynamics [2-4]. One is lead to expect that likewise, radiative 
decays of D mesons [5-10] will bring similar useful insights, as well as providing possible 
checks [11-15] for physics beyond the standard model. 

In a recent paper [|I6|| , we have initiated the study of radiative decays of type D — > 
Kn'-f, with non resonant Ktt. In [[16] we studied the Cabibbo allowed decays which 
have a bremsstrahlung component, D + — > K°7r +r y and D° — > K~7t +/ -f. The theoretical 
framework used is that of the effective Lagrangian which contains both heavy flavor 
symmetry and SU(3)l X SU(3)r chiral symmetry [I7| and the calculation was carried 
out with the factorization approximation for the weak matrix elements. The merging 
of this framework with the known QED requirements for the bremsstrahlung radiation 
necessitated the development of a specific technique for the treatment [|l6j of these decays. 

In the present article we undertake the study of decays with neutral mesons D° — > 
K TT (r],ri')'~f, which are likewise Cabibbo allowed but do not involve a bremsstrahlung 
component. Due to this feature, these purely direct radiative decays may provide a cleaner 
testing ground for the study of the dynamics involved and for exploring the suitability 
of the theoretical framework employed in their calculation. In |l(| it was found that the 
branching ratio for the parity-conserving part of the decay D° — > K~ir +, y, which part is of 
purely "direct" nature, can be as high as 10~ 4 , while the total branching ratio for photon 
energies above 50 MeV is close to 10~ 3 . Such relatively high rates, especially for the direct 
transition, strongly motivates the desirability of studying the decays with neutral mesons 
which we consider here. 

In decays of this type, one does not expect any significant short-distance contribution 
[5-10]. Since the decays D° — > K°P Qr f involve transitions between a rather heavy me- 
son and pseudoscalars, we adopt for our calculation the heavy quark chiral Lagrangian 
(HQxL) containing heavy flavor and SU(3)i x SU(3)r symmetries as the appropriate 
theoretical framework. However, since virtual vector mesons may play an important role 
in these decays in the (Pj) channels, we need to complement the Lagrangian |DJ with the 
light vector mesons. For this purpose we use the generalization of HQ\L by Casalbuoni 
et al. |18 to include the vector mesons in the Lagrangian, in which the original symmetry 
is now broken spontaneously to diagonal SU(3)v (HJ- This framework is described in 
detail in Refs. || |18|1 , henceforth we recapitulate only the main features. The light meson 
part of the Lagrangian is written as 

f 2 

Clight = -'—{HA^ + atr^V^-prf]} 

+ ^tr[F, v (p)F^(p)) . (1) 



with a = 2 for exact vector meson dominance 19 . The vector and axial-vector currents 
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are given by 



V„ = -(vlD^u + uD^) A„ = -ArfDpU - uD^u ] ) , (2) 



where u = exp (jy } , IT being the matrix of the pseudoscalar fields and f = f n = 132 MeV. 
is the covariant derivative. Moreover, p^ = i^P^, with p M being the matrix of the 



vector fields: 



V2 n » \ 



V2 



(3) 



For the vector coupling g v (my) we use experimentally determined values from leptonic 
decays |H| of vector mesons, which accounts for the symmetry breaking effects. 
In the heavy sector, to order 0(p) one has the Lagrangian 

C str = Tr(HJv-D ab H b )+igTr(H a H bllx A^ 5 ), (4) 

where D^ b H b = d^H a — H b V ba , while the trace Tr runs over Dirac indexes. The flavor 
indexes are denoted as a and b. Both the heavy pseudoscalar and the heavy vector meson 
are incorporated in a 4x4 matrix 

H a = ±(l +? f){P„ al »-P 5al5 ) (5) 



The strong coupling constant has been determined [21] recently from D* — > D~n decay to 
be g ~ 0.59 and is the four-velocity of D meson. 

The electromagnetic interaction is introduced by gauging the Lagrangians ([I]) and 
(|j) with the U(l) photon field, thus amending appropriately the covariant derivatives. 
However, the gauging procedure alone does not generate the D* — > D7 transition, thus 
requiring the introduction of additional terms C c . There are two such terms in the frame- 
work used, giving the direct photon-heavy quark interaction with strength A' (being of the 
order 1/mq) and a light vector meson-heavy meson interaction with strength A (being 
of the order 1/A X where A x is the chiral perturbation theory scale). The second term 
provides photon emission via vector meson dominance (VMD). Thus, 

C c = -\'Tr[H a a^F^{B)H a \+i\Tr[H a a, u F^{p) ab H b }. (6) 

An additional contributing term to the radiative decays via VMD is the Wess-Zumino- 
Witten anomalous interaction for the light sector, given by p2| : 



4& = -^^^Tr{d,p v d aP ^). (7) 
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In the calculation, instead of using the SU(3) symmetric coupling, we shall rely on the 



experimentally measured effective couplings of the V — * P7 transitions [23|. For A, A' we 
take A = —0.49 GeV^ 1 , X' = —0.102 GeV~ x , as determined from D* +, ° electromagnetic 
and strong decays JTH], the signs chosen as to conform with the quark model. The effective 



weak Ac = 1 nonleptonic Lagrangian of relevance to the decays investigated here is 

£#/(A C = As = l) = -^V ud V: s [a 1 1 + a 2 2 ], (8) 

where 0\ = {sc) v _ A {ud)^y-A, 2 = {uc) v _ A {sd)^y-A and V^- are the CKM matrix 
elements. The effective Wilson coefficients are taken as a% = 1.26, a 2 = —0.55 J24[]. In our 
calculation we rely on factorization, which implies < P°K°\Ox\D° >= 0. The heavy-light 



weak current is bosonized as 17, 18 



J Q »= § i a Tr[Y{l-l,)H b u\ a ] 

+ a 2 Tr[^ l5 H b v a (p a -V a ) bc u\ a } + ... , (9) 

The constant a is then related to the decay constant fa by a = Jd^^d, while ax, a 2 are 
determined from the values of the form factors appearing in D — ► Vlv decays, as explained 
16 . Their numerical values are then ail = 0.156 GeV 1/2 , \a 2 \ = 0.052 GeV 1/2 . 
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The general Lorentz decomposition of the D° — > K°P 0/ -f amplitude is given by 
Gi 



M = - 



~^VduV: s (F, ((q . e )(p . fc) _ (p . e)(q ■ k)) + F 2 e^e,v a k^) , (10) 



where Fx, F 2 are the electric and magnetic transitions, which are respectively parity- 
violating and parity-conserving. The four-momenta of P°, K and the photon are q, p, k 
respectively and is the polarization vector of the photon. 

In Fig.l we exhibit the Feynman diagrams contributing to the D° — > A'°7r°7 decay. 
The diagrams denoted Ai (parity- violating) derive from the < 7i°\J\D° >< ^°|J|0 > 
term; the same term gives the parity-conserving terms B{, while the parity-conserving 
terms Di come from < 0| J\D° >< 7r°^°| J|0 >. The amplitudes Ai, B>i, Di are related to 
the Fx, F 2 form factors by 

F 1 (D°^K°P° 1 ) = Ax, (11) 
3 2 

F 2 (D° - K°P° 7 ) = J2 B i + J2 D i- (12) 
i=i i=i 

The sums for each row of the amplitudes in Fig.l are gauge invariant. The explicit 
expressions of A i: B iy Di are listed below, for the D° — > K°7r°7 decay. 

„ . v 7 ^ f D f K v ■ k 
Ax = i-^-eg 



2 * / w v-k + v-q + A\vk + A vq + A 
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Figure 1: Feynman diagrams 
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- ief D f K gv^M- 
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(v ■ k)(p ■ k) 

v ■ (q + k) + A) \(q + k) 2 - m 2 + iT p m p ' (g + A;) 2 - m 2 + ir^m, 



(13) 



B 1 = -v^eM^A 



v ■ p 



+ 9: 



v ■ p 



f n \v-k + A (vq + v- k)(v- k + A) (v ■ q + v ■ k)(v ■ q + A) J ' 



+ 9: 
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(14) 



+ 



(q + k) 2 -m 2 + iT p m p (q + k) 2 - m 2 + iY w m u 

M—v-g 
vq+A 



3 V-f. 



—eg K o, K og K *—M — — — ^ — — 

2 Jn {k + P) - m z K * + 11 K *m K * 



2 u-fc + A \ (p + g) 2 — m 2 K * + irriK*^ k* J \ 2 " \3m 



m 



E> 2 = M—e— ( 9p9fmi gg^Wy \ ^ 

2 2 \ (g + k) 2 - m 2 + im p T p (q + k) 2 - m 2 + im w T w ) 



^2 M 3 (m 2 p g ml g^ 
T^JdJk- 



4 (M 2 - m 2 R0 ) \g p (q + k) 2 - m 2 + im p T p g w (q + k) 2 - m 2 + «m w r w / 

The treatment of the decays to r/, 77' requires the inclusion of 77 — 77' mixing. This 
is usually performed [f25| by the mixing of singlet and octet states via a unitary matrix 
characterized by an angle 9. The analysis of various decays involving 77, rf and of meson 
masses has led to a mixing angle in the range (—9°) to (—23°). However, a treatment p6| 
based on chiral perturbation theory and a detailed phenomenological analysis indicates 
that this simple scheme is inadequate and must be extended [Q to include two mixing 
angles 80 and 9g, in addition to the two decay constants. Here we follow this approach with 
1 77 >= cos9 8 \r] 8 > — sin^ol^o > an d \rf >= sin# 8 |ri 8 > +cos6'o|?7o >• For the D° — > K ^ 
and D° — > K i]' / ~f amplitudes the same form as Eqs. ( |13| - |T5D holds, except for replacement 
of constants. In order to obtain these amplitudes one has to replace in above Eqs. f v by 

f v /(lQy/2) or fr//{K2 v/2), while in amplitude D l one replaces ( 1 + , a ) 2 —m 2 K _ 



"K 

,2 



2 
R 

{p+q) 2 —rn 2 K ,+im K *r K * 



by y/2[lQ + (K2 - K' n s ) ( + Z"L m i ) ■ Moreover, one replaces g p7T1 , gw 7 by g pv (r,>)-y and 



\ d 1 V-^d s ''(p+9) 2 -m^ 4 

SW^Ot- The factors Kj, K r J, , are octet-singlet mixing factors, given by 

v _ cos gg sin go / _ sin fl 8 | cos 6 

_ 2 cos 6s _ sin go , ^ 2 sin gg cos6> 

A recent analysis by Feldmann, Kroll and Stech has given further clarification of 
the theoretical basis of this scheme. Following their procedure [f27], PBfl , we use # 8 = —20.2°, 
#0 = —9.2°, = f w and /' = 1.13/^. We have used these values in our calculation, as 
well as checking the sensitivity of the results to the scheme and to the values of above 
constants within an acceptable range. 



Using now Eqs. C|13| - [r5]) we calculate the differential decay distributions and total 
decay rates for D° — ► K 7T 0/ ~f, and with appropriate replacements as indicated above, for 
D° — > K°f](r]')'y. Since p is a wide resonance we use a momentum dependent width 

TM 2 ) = r P «) f ^IV' 2 ^e(4m|) . (18) 



2 -Ami) V¥ 
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We are interested in decays to nonresonant Kir final states and we also wish to delete 
from the final state resonant (K'y) and (Pj) configurations. Accordingly, we subtract 
the calculated rate given by the diagrams that include above configurations from the 
total rate. However, our results include the contributions of the remaining virtual vector 
mesons and their interference with the resonant ones. Our resulting prediction for the 
major radiative decay is 

BR(D° -> K°7r 7 ) NR = 3.0 x KT 4 . (19) 

The parity- conserving and the parity-violating parts of the amplitude contribute approx- 
imately 2/3 and 1/3 to the decay rate respectively. Now, if we do not subtract the direct 
vector mesons contribution, one gets a branching ratio of 3.0 x 10 -3 for this channel. 
Obviously, the two cases have rather different Dalitz plots, the resonances being easily 
identified in it as they dominate the decay rate in the latter case. On the other hand, if we 
exclude vector mesons from our Lagrangian, the total decay rate is reduced to 2.6 x 10 -8 . 
We have also calculated the size of the rate coming from the direct resonant process 
D° — > K* 0/ ~f using our formalism and we find it to be more than two orders of magnitude 
smaller than fll9|), in general agreement with previous estimates || || ||. 

For the decays to 77, 77', again after deleting resonant two body channels from the final 
state, we get: 

BR(D° - KV)nr = 1.1 x lO" 5 , (20) 
BR(D° -> Ky 7 ) NR = 4.3 x 10~ 8 . (21) 

Allowing also 777, 7/7 in the final state, coming from the intermediate D° — > p(oo)j 
decay, these figures are raised to 3.9 x 10~ 5 and 1.4 x 10 -7 respectively. Excluding com- 
pletely the contribution of vector mesons, these branching ratios are lowered to 3.7 x 10~ 8 
and 1.3 x 10~ 9 . In the decays to 77, 77', the dominant contribution to the decay rate is due 
to the parity conserving part of the amplitude. 

The major role played by off-the-mass-shell vector mesons in these decays is evident 
from the above procedure. In particular, the w-meson exchange is the dominant contri- 
bution to the pionic decay (Eq.(20)), while for the decay D° — > K i]'~f, both u and K* 
exchange are giving major contributions. For the decay to 7/, the main contribution is 
due to K* . The contribution of the p meson is generally much smaller, in all channels. 

We have also found that the D° — > Krj'j rate is smaller than D° — > K^^j rate partly 
due to phase space, but mostly due to the different values of the coupling constants 
involved. The smallness of the decay with rj in the final state is mainly due to phase 
space inhibition. 

In Fig. 2 we display the Dalitz plots for the three decay channels, and the photon 
spectra of these decays. The Dalitz plot coordinates are m 12 = J{P — k) 2 on x axis and 

m 23 — y{P~~ p) 2 011 y ax is an d the spectrum is expressed in terms of (m 12 ) 2 = (M — k) 2 . 
In the main channel D° —>■ K TT°~f the peak of the parity-conserving contribution is at 
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E 1 = 700 MeV and of the parity-violating contribution at E 1 = 500 MeV. The total 
spectrum is predicted to peak around E 1 = 650 MeV. For the rj decay, the spectrum 
peaks at £ 7 = 300 MeV. 

Before concluding we wish to make some observations concerning the appropriateness 
of the theoretical framework we used. By using the Lagrangians (|), ©-(@) as well as 
the factorization approximation we can calculate also various decays of type D — > VP. 
In particular, we have calculated the rates for four channels of this type, which relate 
to the radiative decays we considered. Of these, u channel is the most important one 
for the pionic decay while for the rj channel, both the u and K* channels are giving the 
main contribution. We present now the results we obtained, giving in parentheses the 
observed branching ratios: BR(D° -> cuK°) = 5.6%(2.1%), BR(D° -> pK°) = 7%(1.2%), 
BR(D° -> ttK*) = 11%(3.2%), BR(D° -> r]K*) = 3.2%(1.9%). The rates we obtain are 
within factor two-four (the relevant ones are on the lower side of this) of the experimental 
ones, which leads us to estimate that this is essentially the accuracy of our calculation, 
including also the uncertainty due to subtraction of the resonant channels from the total 
rate. We have checked the sensitivity of our results to the values of the mixing angles and 
the values of f v , f . Within reasonable values for these constants, the changes in rates 
are negligible. 

We conclude by pointing out that the branching ratio of D° — > K°7r°'y is quite large, 
due to the contribution of light virtual vector mesons. That makes it rather appealing for 
the experimental study since in this decay mode there is no need for special treatment of 
the bremsstrahlung component. 

With the abundance of D's produced at B factories, at Tevatron and at the forthcoming 
charm factories, we look forward to experimental results for these decays, especially for 
the 7r and rj modes, which are predicted to have sizable rates of 3 x 10~ 4 and 1.1 x 10~ 5 
respectively. 
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